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We study fermion reflection at a phase wave which is formed during a 
bubble collision in a first order phase transition. We calculate the reflection and 
the transmission coefficients by solving the Dirac equation with the phase wave 
background. Using the results we analyze the damping and the velocity of the 
wave. 
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There have been many studies on the cosmological roles of the first-order 
phase transition which proceeds by nucleation and collisions of vacuum bubbles 0, 
especially in some inflation models 0. 

If a global or a local symmetry is broken at the transition, there could be a 
phase wave, which has not been studied fully. In this paper, we study a fermion 
scattering at the phase wave. 

At the nucleation, due to a finite correlation length, the two true vacuum 
bubbles that are separated by more than the correlation length or the horizon 
distance may have different phases. Such a possibility is a key ingredient for the 
Kibble mechanism f3|]. 

If there is a phase difference between the two colliding bubbles, a pair of 
phase waves occur and propagate into each bubble. The phase wave [|, || is a 
mechanism by which the bubbles get their new phase values determined by the 
so-called "geodesic rule" in a first || or a second Q order phase transition. 

The energy of the colliding bubble wall turns into that of the modulus wall 
and the phase wave. If the phase difference of the two bubbles is of order 1, the 
phase wave carries away most of the energy of the bubble wall ||. So the phase 
wave can be a very energetic phenomenon as the bubble wall is. 

Furthermore, investigating the particle scattering at a moving bubble wall is 
important to study the bubble kinetics and cosmology. For example, to calculate 
the velocity and the width of the electro- weak bubbles [0] and the CP violat- 
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ing charge transport rate by the wall for baryogenesis ||, one should know the 
reflection coefficient of fermions (e.g. top quarks) scattering at the wall. 

However, there have been a few works about the dynamics and cosmological 
roles of the phase wave compared with the bubble walls. 

So, in this paper, we investigate the interaction between the phase wave and 
fermions. In the cosmological phase transition, a propagating phase wave may 
collide with the fermions in the plasma. 

Consider a complex scalar field 4>(x) = p(x)exp(i9(x)) whose Lagrangian den- 
sity is given by 

C= 1 -d,<j>*d^-V(p), (1) 

where V(p) is a potential appropriate for a first-order phase transition. 

Since the potential depends only on the amplitude p, one may choose an 
arbitrary phase for the minimum of the potential, (0). Of course, this is due to 
the global U(l) symmetry the theory possesses. In general, we can make the true 
vacuum to lie along the direction of the real part of 0, because the phase of <fi 
can be absorbed by redefining the phase of the fermion field. However, if (</>) is 
space-time dependent, we can not globally rotate out the phase which may play 
physical roles. 

From £ we obtain the equation of motion for 9; 

d^df.9 + -^9d^p = 0. (2) 

Without loss of generality, we can choose the two phases for the colliding 
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bubbles as zero and 2A6*(< n), respectively . Then, by the geodesic rule, the 
overlapping region of the two bubbles should have A6(< tt/2). 

We consider the case where the width of the phase wave (d) is much smaller 
than the radius of curvature (r) of the phase wave. In this case we can treat our 
problem as one dimensional. We choose z axis as the direction normal to the 
plane of the phase wave, and assume that momentum of the incident fermion is 
parallel to the z axis. The solution for the more general momentum can be found 
by appropriately Lorentz boosting the one dimensional solution. 

A possible solution of eq.(Q) is 

9(z, t) = A0[1 + l -f{z + vt) + l -f{z - vt)], (3) 

where f(z) is some profile function which has an asymptotic value ±1 at z = ±00, 
respectively. 

The second term indicates the left moving phase wave(— z direction) with 
velocity v relative to the plasma, while the third term is the right moving one 
(+z direction). Since d^p = for the phase wave, v is the speed of light c. 
However, eq.(Q) ignored the pressure by the fermion scattering which is just what 
we are going to study in this paper. So we assume \v \ < c. 

To calculate the reflection coefficient of the fermion we choose the rest frame 
of the left moving phase wave (v = 0) and assume that the incident fermion 
comes from the left. (See Fig.l.) When the typical wave length and the mean 
free path of the fermion are much larger than the phase wave width ~ 1/wio, we 
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can approximate f(z) as 2Q(z) — 1 with step function 0(z),i.e., 9{z) = A8 when 
z > and zero for z < 0. 

Consider a Lagrangian density for the fermion \1/ 

= + ^! r ^r - (h^ L ^ R (f) + h*ty R ^ L (f)*) (4) 

with the (real) Yukawa coupling constant h. In the rest frame of the phase wave, 
^ acquires position dependent mass m(z) = h < <p(z) >. Therefore we get the 
equation of motion for \1/ from L^; 

(tf-m(x)P R -m*(x)P L )V = 0, (5) 

where Pr^l are the chirality projection operators. 
We adopt the following ansatz||, 



V(t,z) = (i$ + m*(z)P R + m(z)P L )e- iaEt ij(k,z) 

= (a£ 7 ° + i^d z + m R (z) - tmj(z) l5 )e- tC7Et ^(k } z), (6) 

where a field ip satisfying following the Klein-Gordon-like equation is intro- 
duced; 

(E 2 + d 2 z - \m(z)\ 2 + im' R (z)j 3 - m'^z)^)^) = 0. (7) 

Here prime denotes the derivative with z, m(z) = m R (z) + imj(z) and a is +1 
and —1 for fermion and antifermion, respectively. 

If m'(z) ^ 0, there is an additional potential term proportional to the spatial 
change of m(z), which is the origin of the scattering. For our case, there is a 
delta function type potential at z = 0. 



Let us find momentum eigen states. We expand ip in the eigenspinors of 
7 3 with eigenvalue ±z satisfying 7 3 Mj_ = ±iu±, where u\. = 1/Vz(l, 0, ±z, 0) T 
and w± = l/\/2(0, 1, 0, =F«) T with spin index s = 1,2. (We use here the same 
conventions for the 7 matrices and u s ± in ref.||.) Note that 7°u^_ = 

Since m(z) = moe ie ^ z \ ttir = mo and mj = for z < (region J) and 
mn(z) = mocos(A6) = m c and mi(z) = mosin(A8) = m s for z > (region //). 

From now on we will consider the case of particle(cr = 1). For z 7^ 0, m'(z) = 
and eq.(0) becomes a free Klein-Gordon equation. Then, the most general right 
moving solution to eq.(^) can be represented as a combination of u s + and «i,.i.e., 

^ = ^£53 <?•«•, (8) 

s ± 

where k = ^ E 2 — m\ and C± are some constants. 
Inserting if) into eq.(||) we get \I/ for 3 < 0; 

$ = g-ifit+ife^ £{C£(£i£ + (mo - iA:)«;) + Ci(Ew; + (mo + iJfe)«L)}- (9) 

s 

These terms are not independent in each other in u± basis, so we can set Ct = 
without loss of generality. Therefore, 

ifj(z)=e lkz Y.C> S + - (10) 

s 

Similarly, for the left moving wave we can choose \I/ proportional to (Eu s + + (m — 
ik)ui_). So for the region /, the general solution can be given by 

= e~ iEt C s + {e lkz [{m - ik)u s + + Eu s _J\ 

s 

+R s e~ ikz [Eu s + + (m - ik)u s _]}. (11) 
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The first term is a incoming wave and the second term is a reflected one. In the 
region II only a transmitted wave exists; 

* w = e~ lEt+lkz £ T S C;{(£ - m s (-Y)u s _ + (m c - (12) 

Here, i? s and T s are some spin-dependent coefficients. 

Now, equating and at z = and comparing each coefficients of we 
get two algebraic equations for R s and T s ; 

-E + (E - (-) s m s )T s 



Rs 



rriQ — ik 



m mn — ik + ER S , , 

T s = — r, ~, 13 

m c — ik 



which has a solution 



R s = — [E(m -m c )-(-) s m s (m -ik)], 

T s = ^[(m -ik) 2 -E 2 }, (14) 

where D = (m c — ik)(m Q — ik) — EiE — (— ) s m s ). 

From and ^>u we get a vector current j 3 = ^7 3 \]/ for the region / and //. 

jf = 2kEj2\C-l\ 2 (l-\R s \ 2 ), 

s 

f n = 2kEY,\Cl\ 2 \T s ni-t)^±y ( i5) 

E 

During this calculation the orthogonality condition u s ± u s ^ = 5 SS > and u±u± = 

are useful. Here, from eq.(14) 

IP ,2 2m (m -m c )(l-(-)'^) 
\R 8 \ - j , 

\T S \ 2 = ~y (16) 



with F = 2(2k 2 + mo(mo — m c ))(l — (—) s m s /E). We have used the relations 
E 2 = k 2 + wig and m 2 + m 2 , = m\ repeatedly. 

So the reflection and the transmission coefficients are given by 



R 



T 



m (m - m c 



^2\C+\ 2 2k 2 + m (mo — m c ) ' 



^|C^| 2 2k 2 + m (m - m c ) 



(17) 



It is clear that the unitarity follows from eq. (p77|) , i.e., R + T = 1. We can get 
the coefficient for an antifermion by substituting — E, —k and (corresponding 
to C s + ) instead of E, k and C+ respectively. Since R and T are not changed by 
these substitutions, we can say that the fermion and the antifermion reflect at the 
phase wave in a CP conserving way. It is contrary to the bubble wall case where 
a position dependent phase of generally induces a CP violating reflection. 

Since R and T depend on cos(A8), there is no difference between A# and 
— AO. It also means that R and T are not dependent on the incoming direction 
of the fermion (left or right). 

Now let us find approximate R for A# <C 1. In this case, m c ~ mo(l — A6 2 /2) 
and m s ~ m A9. Gathering all terms up to 0(A9 2 ) we find 

mlAQ 2 , , 



R 



4A; 2 v ' 

This implies that the particles with momentum k k c = m§AQj2 will totally 



reflect at the phase wave and transfer the momentum 2k to the wave. 

We also confirm this result fllSD by solving eq. (H) and eq. (|7|) in a second-order 



perturbation [1 | when the width of the wave d is finite. There we checked that R 
in eq.(|l8D is the leading order term at d —>■ limit, expanding e %e ~ 1 + iO — 8 2 /2 
and finding solution up to the second order, 

From eq.(|T?D we find that the definition for k c is also reasonable even for 
A0 ~ 1. Hence we will approximate R(k) as 8(k c — k). 

When there is no viscosity, as the phase wave expands, it should lose its 
gradient energy {d^O) 2 and its width d becomes thicker and thicker, because, 
contrary to the bubble wall case, there is no net pressure from the vacuum energy 
difference across the phase wave[|J. 

So it is possible that the particle scattering may reduce the velocity of the 
wave and even stop it. 

The viscosity by this reflection could be calculated by the thermally averaged 
momentum transfer (P). 

From the energy momentum relation E 2 = k\-\-k 2 -\-m^ where k\ = k 2 +k 2 , we 
find k z dk z = EdE. In thermal equilibrium states, there are n(E)d 3 k/(27r) 3 \k z /E\ 



particles per unit wall area in unit time providing a momentum transfer 2k z W2 



Here n(E) = [e E l T + 1) 1 is the Fermi distribution function at the temperature 
T. 
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Assuming v > 0, we obtain 



/rp' I? . /"OO 
dfc,2fc,i2(fc,)-^n[ 7 (^ + 

+ /° ^2£^(£g^7i[ 7 (£ + *;£;,)]}. (19) 

J-OO III 

where d 2 kj_ = dk x dk y . The first term corresponds to the pressure by the fermions 
with k z > 0, while the second term corresponds to that with k z < 0. Changing 
integration variable from k z to — k z in the second term we obtain 

/d^k i r°° k 2 
J o dk z 26(k c - k z )^{nYi{E + vk z )} - n [y(E - vk z )]} (20) 

Fig. 2. displays a numerical result of this integration. (P) increases linearly as v 
increases, and then goes to zero as v approaches to 1. 

Now we assume that the velocity of the wave is nonrelativistic (v <C 1). Then 
n(rf(E + vk z )) ~ n(E) + vk z ^§^ = n(E) + v(3k z n(E)(n(E) - 1). Changing 
integration variables from d 3 k to 2irEdEdk z and assuming (3E 3> 1 give 

P [°° Jrn _ 0E fV^< 



H roc f\ & m o 

(P) = -«4 / dEe~ pE / dk z 9{k c - k z )kl 

H Jmo JO 



ft fE c r JE 2 -ml 

P { / dEe~ pE / dk z 9(k c - k z )k\ 

Jmo JO 
roc rk c rW E 2 —ml 

\ dEe' pE [ dk z 6(k c -k z )kl+ dk z 6(k c -k z )k 3 z ]}, 
Je c Jo Jk c 



= ~ v 2 - . 



+ 



(21) 



where E 1 ^ = + k\ — ml(l + A# 2 /4) and we used the approximation n(E) 
exp(—(3E). We have separated the integration regions to drop the third term. 
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Then, 

TP 



(P) ~ - V JL f c d Ee-^(E 2 -m 2 ) 2 -v^-e'^ 

4:7l 2 Jm 47T 2 

e~ pE % (22) 



471-2 

because the first term is of order (A#) 6 , while the second term is of order (A#) 4 . 
For 9 = 9{r — vt) the Newtonian equation for r is given by[]I3] 

(Pr 2a , , 

= ~ + <P> (23) 

where v is now dr/dt. This equation can be derived by considering the force acting 
on the unit spherical area of the wave with the 'radius' r and 'mass' density a. 
The surface energy density a is related to the energy momentum tensor by 

a = JdzT° = /4(^) 2 ^. (24) 
Hence for nonrelativistic and planar approximation^ — > oo) 

which has a solution 

t> = Voexp[—- — ^—exp(—(3E c )t}, (26) 

where vq is an initial velocity. 

Hence the phase wave stops in a time scale t stop = ^ E gr^ iEc and the wave 
travels the distance 

Ar ~ v t stop [l - exp{-t/ tstop)} (27) 
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during t. 

The bubble nucleation rate per unit time and volume is T ~ T 4 e~ A and 
A ~ 4inMp/T0. The average bubble distance and the bubble percolation time 
scale is t bub ~ where j3 = ^AH. Here the potential form dependent factor 
7 = dlnA/dlnT is generally 0(1) ~ O(10 2 ) p[. 

For example, if T ~ 10 2 GeV, A ~ 10 2 and t buh ~ lO 2 ^ 1 ~ 10 2 M P /T 2 . On 
the other hand, if all the coupling constants are not too smaller than unity, then 
m ~ T, and t stop ~ T~ l e f3Ec . Hence, for (3 not too smaller than E c , t stop should 
be shorter than tbub- Generally if T e~ l3Ec Mp ) the phase wave may stop before 
it reaches the other side of the bubble wall. 

If these stopped phase waves persist long time, their energy density may give 
rise to a problem similar to the domain wall problem 0. Even if some ignored 
mechanisms ( for example surface tension) make the wave disappear, they take 
a role only within the horizon at that time and hence the problem could not 
removed completely So in this view point, it might be preferable to have a 
non-degenerate unique vacuum for (0). 

When the Yukawa coupling h is complex, defining m c = mo(Re(h)cos(A6) — 
Im(h)sin(A6)) and m s = m Q (Re(h)sin(A8) + Im(h)cos(A6)) gives the same 
result in eq.(|T^). 

If there is another scalar field x interacting with <j) through (A/2)A|0| 2 |x| 2 , X 
gets a constant mass \ l ^ 2 p and satisfies Klein-Gordon equation. Hence the scalar 
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field x does not reflect at the phase wave. 

The origin of fermion scattering at the phase wave is a position dependent 
phase change of its mass, while that of the bubble wall is a change of absolute 
mass and, therefore, a momentum change across the wall. 

Our results may be also applied to detection of a spatial phase change and its 
velocity v, if any, in the universe or a condensed matter. For example, by using 
eq.(|T8|) and observing the difference between the two i?'s for electron beams with 
the opposite incoming directions, (say, from region I and II in fig.l., respectively.) 
Each electron beams should have different k in the phase wave rest frame moving 
with v , and have different R. 

The results, so far, are for a global symmetry without gauge fields. Including 
gauge fields may give more complicated and interesting results. 

In summary, we calculated the reflection coefficient of the fermion scattering 
off the phase wave, which could be a first step toward a comprehensive under- 
standing of the phase wave dynamics in the hot plasma. As an example, we 
investigated the damping and the velocity of the wave due to this scattering, and 
found that this damping could even stop the wave in some cases. 

One of the authors (Lee) is thankful to H.M. Kwon, J.D. Kim, Y.D. Han and 
H.C. Kim for helpful discussions. This work was supported in part by KOSEF. 
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FIG.l. 

9(z) in the phase wave rest frame. \&j consists of the incoming and the reflected 
waves and ^ n is the transmitted wave. 

FIG.2. 

The damping pressure (P) in units of Wq versus the velocity of the phase wave v 
for A6 = 0.1 and T = 2mo- The solid curve shows the results of eq.(PD|) and the 
dotted line shows the results of eq. (|22|) . 
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